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pL( " Abstract 

We study a noninteracting supersymmetric model in de Sitter spacetime. A soft 
supersymmetry breaking induces a nonzero vacuum energy density. A short distance 
cut-off of the order of Planck length provides a matching between the vacuum energy 
density and the cosmological constant related to the de Sitter expansion parameter. 
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It is generally accepted that the cosmological constant term which was introduced ad-hoc 



in the Einstein-Hilbert action is actually related to the vacuum energy density of matter 
fields. The vacuum energy density estimated in a simple quantum field theory is by about 
120 orders of magnitude larger than the value required by astrophysical and cosmological 
observations. On the other hand, in a field theory with exact supersymmetry the vacuum 
energy, and hence the cosmological constant (CC), is equal to zero as the contributions of 
fermions and bosons to the vacuum energy precisely cancel pQ. 

A nonzero CC implies the de Sitter symmetry group of spacetime rather than the Poincare 
group which is the spacetime symmetry group of an exact supersymmetry. Based on observa- 
tional evidence for an accelerating expansion [2j[3lll], the vacuum energy density dominates 
the total energy density today. Hence, the spacetime today is close to de Sitter approaching 
asymptotically a de Sitter universe with metric 

ds 2 = dt 2 - a{t) 2 dx 2 } (1) 

where a = e Ht and if is a constant. This metric describes empty space with cosmological 
constant A = 87iGp\, where the vacuum energy density p\ is related to H by the Friedman 
equation: 

H* = ( "V = M pA . (2) 



mJ 3 

The structure of de Sitter spacetime automatically breaks the supersymmetry [5]. Con- 
versely, a low energy supersymmetry breaking could in principle generate a nonzero CC 
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of an acceptable magnitude. Unfortunately, the scale of supersymmetry breaking required 
by the particle physics phenomenology must be of the order of 1 TeV or larger implying a 
CC too large by about 60 orders of magnitude. However, some nonsupersymmetric models 
with equal number of boson and fermion degrees of freedom have been constructed [6] so 
that all the divergent contributions to the vacuum energy density cancel and a small finite 
contribution can be made comparable with the observed value of the CC. 

In this paper we investigate the fate of vacuum energy when an unbroken supersymmetric 
model is embedded in asymptotically de Sitter spacetime. We do not claim that our model 
describes a realistic scenario but it is tempting to speculate along the the following lines. 
Large scale observations reveal that dark and baryonic matter gravitationally cluster occu- 
pying relatively small volume of space in the form of poor and reach clusters connected by 
filaments and sheets [TJ. Most of the volume is occupied by large scale voids almost empty 
of both baryon and dark matter. 

Our working assumption is that voids contain no matter apart from fluctuations of a 
supersymmetric vacuum as a relict of symmetry breaking in the early universe. The early 
universe with exact supersymmetry underwent a set of symmetry breaking phase transi- 
tions. Before the supersymmetry breaking, domains with different vacua may have been 
formed with domain walls separating the domains [Sj. Then, at a later time, supersymmetry 
breaking took place in some if the domains, which thereafter remained populated by dark 
and baryonic matter. Our main assumption is that supersymmetry remains unbroken in the 
voids. However, since the global geometry is de Sitter, the lack of Poincare symmetry will lift 
the Fermi-Bose degeneracy and the energy density of vacuum fluctuations will be nonzero. 
This type of "soft" supersymmetry breaking is similar to the supersymmetry breaking at 
finite temperature where the Fermi-Bose degeneracy is lifted by quantum statistics [9j [10] . 

The remainder of the paper is organized as follows. In section [5] we describe our model. 
The calculations and results are presented in section [3] and concluding remarks in section HJ 



2 The model 

Here we consider a noninteracting Wess-Zumino supersymmetric model with N species and 
calculate the energy density of vacuum fluctuations in de Sitter spacetime. In general, the 
supersymmetric Lagrangian C for N chiral superfields has the form [11] 

£ = X>t$.| D + w($)| F + h.c., (3) 

i 

where the index % distinguishes the various left chiral superfields $j and H / ( ( I ) ) denotes the 
superpotential for which we take 

W(*) = ^ J Em i *&i. (4) 

i 

Eliminating auxiliary fields by equations of motion the Lagrangian (j3j may be recast in the 
form 

C = d^jd^i - m^| 2 + ^a^i - i% , (5) 
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where 0, are the complex scalar and the Majorana spinor fields. For simplicity, from now 
on we suppress the dependence on the species index %. 

Next we assume a curved background spacetime geometry with metric g^. Spinors in 
curved spacetime are conveniently treated using the so called vierbein formalism. The metric 
is decomposed as 

9v»{x) = e a ll e b v ri ab , (6) 
where the set of coefficients e a M is called the vierbein. The action may be written as 

S = J d 4 x^g-(£ B + C F ) } (7) 

where Lb and Lp are the boson and fermion Lagrangians, respectively. The Lagrangian for 
a complex scalar field may be expressed as the sum of the Lagrangians for two real fields 

^ = iE(^v>:,-^v 2 )- (8) 
z i=i 

The fermion part is given by [13] 



C F = % - (#7^ - *;//r^) " \™$rtt, (9) 
where 7 M are the curved spacetime gamma matrices 

r = (io) 

with ordinary Dirac gamma matrices denoted by j a , and e a M is the inverse of the vierbein. 
The covariant derivatives of the spinor are defined as 

* ;M = * rf ,-r M * J (n) 



where 



n • ^r,,- (12) 
r, = ^/V,7 6 ], (13) 



with the spin connection 

< 6 = ^V(e%-r;/ A ). (14) 

In FRW metric the vierbein is diagonal and in spatially flat FRW spacetime takes a simple 
form 

e a M = diag(l, a, a, a) (15) 
where a = a(t) is the cosmological expansion scale. 
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3 Calculation of the vacuum energy density 

It is convenient to work in the conformal frame with metric 

ds 2 = a(r]f(drj 2 -dx 2 ). (16) 
where the proper time t of the isotropic observers is related to the conformal time r\ as 

dt = a(rj)drj. (17) 
In particular, we will be interested in de Sitter spacetime with 

a = e«' = -±. (18) 

In order to calculate the energy density of the vacuum fluctuations we need the vacuum 
expectation value of the Hamiltonian. The Hamiltonian may be expressed as the sum of the 
boson and fermion parts 

n = n B + n F . (19) 

From (EH9]) with metric (1161) we obtain 

He = E ( V) 2 + a^W) 2 + ™V 2 ) > ( 20 ) 

U F = -i-^ {fi'djij} - (9^)7^) + 2~3 m ^- ( 21 ) 

Consider first the contribution of the scalar fields. Each real scalar field operator is 
decomposed as 

V ( V ,x) = ^a" 1 ( X k(v)e rk "a k + X k(vY ^ ? 4) , (22) 

fe 

where a k and a|. are the annihilation and creation operators, respectively. The function \k 
satisfies the field equation 

X'l + (m 2 a 2 + k 2 - a"/a) Xk = 0, (23) 

where ' denotes a derivative with respect to the conformal time rj. In massless case the exact 
solutions to this equation may easily be found [13]. In particular, in de Sitter spacetime 
a" I a = 1/r] 2 , and one finds positive frequency solutions 



* = vk e ~' tv (' " i) • (24) 

The operators a k associated to these solutions annihilate the adiabatic vacuum in the asymp- 
totic past (Bunch-Davies vacuum) [T3| fT5]. 

If m 7^ solutions to ( 1231 may be constructed by making use of the WKB ansatz 

Xk (r,) = 1 e-T'^W* (25) 
/2VaW k (v) 



where the function Wk may be found by solving (1231) iteratively up to an arbitrary order in 
adiabatic expansion [12]. For our purpose we need the solution up to the 2nd order only 
which reads 



W k = oj k + u 



(2) 



where 



Uk = Jm 2 + k 2 /a 2 . 
The general expression for the second order term is [12] 



U) 



(2) 



3 1a 2 3 1a 3 k 2 a 2 1 k 2 a 5 k 4 a 2 



> Wfc a^ 



4 w fc a 4 a 2 ul a 2 ' 4 a 2 ^ a ' 8 a 4 a;| a 2 



(26) 
(27) 

(28) 



where overdot denotes a derivative with respect to t. For de Sitter spacetime, using ( ITS"]) we 
obtain 

# 2 



1 + 0(m 2 /u 



(29) 



We can calculate now the vacuum expectation value of the boson Hamiltonian. Using 
the properties of a and a) and replacing the sum over momenta by an integral in the usual 

way 



from (1221) we find 

>= 

Using 023]) with §M) we obtain 



7TT 



(2tt)s 



<U B > 



d 3 k 



u k + 



IIP 

2 u k 



1 + 



m 



H 2 



(31) 



(32) 



The first term in square brackets is identical to the flat spacetime result. The second term 
is a quadratically divergent contribution due to de Sitter geometry, the next two terms are 
logarithmically divergent, and the rest is finite. 

Next we proceed to quantize the fermions. The Dirac equation in curved spacetime may 
be derived from (Q. Specifically for a spatially flat FRW metric we obtain 



3 a 



2a 



il° (d + — * + i-y'a* - mtf = 0. 



(33) 



It is convenient to rescale the Majorana fermion field \1/ as 

= a -' i/2 i). (34) 
and introducing the conformal time we obtain for ip the usual flat spacetime Dirac equation 



i^drjip + irfdjij) — amip = 0, 



(35) 
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with time dependent effective mass am. The quantization of ip is now straightforward [TBI [IT]. 
The Majorana field ip may be decomposed as usual 

1>(ri, *) = E {u ks (v)e rk %s + v ks (v)e- rk % s ) , (36) 

k,s 

where the spinor Uk s may be expressed as 

1 ( {Kit + am(k)<Ps \ 

nt - = 7v{ 3h„.4>, )■ (37) 

Here, the two-spinors <p s are the helicity eigenstates which may be chosen as 



Ih *-=(! 



(38) 



The spinor Vk s is related to by charge conjugation 

Vks = ^V^f- (39) 

The norm of the spinors may be easily calculated 

I i 

UksUks = -VksVks = y(amQ - <if)( am Cfc + K'k) ~ y k<2 \Ck\ 2 - (40) 

The mode functions (k satisfy the equation 

Ck + (m 2 a 2 + k 2 - ima')Ck = 0. (41) 
In addition, the functions (k satisfy the condition [17] 

k 2 \( k \ 2 + (amQ - Kk)(amCk + <l) = C\. (42) 

It may be easily verified that the left-hand side of this equation is a constant of motion of 
equation (|4ip . The constant C\ is fixed by the normalization of the spinors and by the initial 
conditions. A natural assumption is that at t — (77 = — 1/H, a — 1) the solution behaves 
as a plane wave ( k = C 2 e~ lE \ where E = y/k 2 + m 2 . This gives Cfc(0) = C 2 , Ck(0) = —iC 2 E, 
and hence C\ = 2C 2 E(m + E). From and (325 we obtain 

UksUks = -VksVks = ^(C 2 - 2A; 2 |C fc | 2 ), (43) 

which at t = reads 

777; 

MfcsWfcs = -V ks Vks = C lyg ■ ( 44 ) 

For C\ = 1 this coincides with the standard flat spacetime normalization [T3] . 

In massless case the solutions to (14ip are plane waves. For m ^ two methods have 
been used to solve ( ]4il for a general spatially flat FRW spacetime: a) expanding in negative 
powers of y/ m 2 + k 2 and solving a recursive set of differential equations [16] b) using a WKB 
ansatz similar to ( )25|) and the adiabatic expansion |17j . 
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By making use of the decomposition (!36|) and the standard anti-commuting properties 
of the creation and annihilation operators, the vacuum expectation value of the fermion 
Hamiltonian (J2T1) may be written as 



< H F >= Yl 53 v >°s ( am ~ k ^)v ks . 

k,s 



(45) 



Evaluating the expression under the sum and replacing the sum with an integral as in (130]) 
we obtain 



< U F >= -7 



d 3 k 



(46) 



The expression in square brackets was calculated in [16] for a spatially flat FRW metric. We 
quote their result for the divergent contribution: 



< Hf >div = 



d 3 k 



—E 



l)m 2 



2E 



+ 



l) 2 m 4 



8E 3 



+ 



{a'fm 2 
8E 3 



(47) 



Note that the first three terms in square brackets are identical to the first three terms in the 

m 2 in powers of E~ 2 . Hence we can write 



expansion of auj k = \JE 2 + 



< T-Lf >div— ~7 



d 3 k 
(2^p 



{a'fm 2 
8a*u 3 



(48) 



The first term in square brackets is precisely the flat spacetime vacuum energy of the fermion 
field. The second term is a logarithmically divergent contribution due to the FRW geometry 
and the last term is finite and vanishes in the flat-spacetime limit a' — > 0. Note that, as 
opposed to bosons, there is no quadratic divergence of the type H 2 /ujk- 

Assembling the boson and fermion contributions, the final expression for the vacuum 
energy density of each chiral supermultiplet is 



Pv 



-< Hf + 7~Lb >— 



d 3 k 1H 2 
(2^)32^ 



5 m 



H 2 



(49) 



The dominant contribution comes from the first term which diverges quadratically. To make 
the result finite we change the integration variable to the physical momentum p = k/a and 
introduce a cutoff of the order of the Planck mass A cut ~ mpi. The leading term yields the 
total energy density of the vacuum fluctuations 

NH 2 /-Acut , _ r . NH 2 K 2 



Pv 



4tt 2 



V 



dp (l + C{p- 



cut 



8tt 2 



(l + 0(A~ C 



In A 



cut 



(50) 



where N is the number of chiral species. Note that the leading term in ( 149~|) is due to 
bosons; fermions only provide a precise cancellation of the quartically divergent flat spacetime 
vacuum term. 

If the Einstein cosmological term is precisely zero then the only source of the vacuum 
energy density are the vacuum fluctuations of matter fields. Hence we identify p vac = p\ 
and if we compare (J50l) with the Friedman equation (j2]) we find that our cutoff should satisfy 



A,. 




mpi . 



(51) 
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It is worthwhile to note that several approaches [IH1 UHl [20l [2H [22l [23l [231 HH] with 
substantially different underlying philosophy lead to results similar to (!50|) . Fluctuations in 
the energy density of a scalar field were estimated [21] yielding the same order of magnitude 
as in ( 1501) . There, a cosmological horizon radius Rh = 1/H was employed as a long distance 
cutoff. This idea is similar in spirit with [22] where an upper bound 

Pvac^AfW<|^ (52) 

was proposed from a holographic principle. Here, Auv and L denote the ultraviolet and 
long distance cutoffs, respectively. Our result would saturate the holographic bound (152]) 
provided L = 1/H. 

In recent papers (23J, E3] a residual quadratic contribution of the form H 2 Al ut has been 
found after canceling the flat spacetime contribution. In particular, the work [23] presents a 
calculation of the zero-point energy using only a massless boson field and obtains two types 
of contributions: the quadratic and the quartic type A^ ut . Then, the quartic contributions 
to CC was simply canceled by hand on the basis of the procedure used previously in the 
literature with the so-called ADM mass. As we have demonstrated here (see also [25]), in a 
supersymmetric world such a cancellation by fiat is unnecessary because the cancellation be- 
tween bosons and fermions of all (not only quartically divergent) flat-spacetime contributions 
is naturally provided by supersymmetry. 

Another important point [23l |25] is that the vacuum fluctuations cannot be interpreted 
as a part of CC because the vacuum fluctuations do not yield the equation of state p = —p, 
as a consequence of the energy momentum tensor not having a CC form. This behavior 
was already observed in flat space time if a three dimensional cutoff regularization was 
employed [26]. A covariant regularization in flat space time should yield the vacuum energy 
momentum tensor of the form = pf]^ u - Naively, in curved spacetime one would generalize 
this to the CC form = pg^ v . However, since a curved geometry involves the Riemann 
tensor and its covariant derivatives we may expect the energy momentum tensor at linear 
curvature order to be of the form T^ u = (a + (3R)g flu + iR^u where a, (3, and 7 are constants 
that do not depend on curvature. In reference [27] the vacuum contribution to the energy 
momentum tensor has been investigated using an explicitly covariant regularization scheme 
for an unbroken supersymmetric model embedded in a general curved spacetime. One loop 
contributions to the effective potential were calculated using a covariant UV cutoff in an 
approach similar to Sobreira et al [28J. 

4 Conclusion 

We have found that the leading term in the energy density of vacuum fluctuations is of 
the order H 2 vn^ x if we impose a short distance cutoff of the order . In this way, if we 
require that the de Sitter expansion parameter H equals the Hubble parameter today, the 
model provides a phenomenologically acceptable value of the vacuum energy. We have also 
found that a consistency with the Friedman equation implies that a natural cutoff must be 
inversely proportional to y/N. A similar natural cutoff has been recently proposed in order 
to resolve the so called species problem of black- hole entropy [29] . 
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